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Abstract 

In this paper, I address the oscillation probability of O(GeV) neu- 
trinos of all active flavours produced inside the Sun and detected at the 
Earth. Flavours other than electron-type neutrinos may be produced, 
for example, by the annihilation of WIMPs which may be trapped in- 
side the Sun. In the GeV energy regime, matter effects are important 
both for the "1-3" system and the "1-2" system, and for different 
neutrino mass hierarchies. A numerical scan of the multidimensional 
three-flavour parameter space is performed, "inspired" by the current 
experimental situation. One important result is that, in the three- 
flavour oscillation case, P a p ^ Pp a for a significant portion of the 
parameter space, even if there is no CP-violating phase in the MNS 
matrix. Furthermore, P^ has a significantly different behaviour from 
P TT , which may affect expectations for the number of events detected 
at large neutrino telescopes. 



1 Introduction 



In the Standard Model of particle physics, neutrinos are strictly massless. 
Any evidence for neutrino masses would, therefore, imply physics beyond 
the Standard Model. Even though the direct experimental determination 
of a neutrino mass is (probably) far beyond the current experimental reach, 
experiments have been able to obtain indirect, and recently very strong, 
evidence for neutrino masses, via neutrino oscillations. 

The key evidence for neutrino oscillations comes from the angular depen- 
dent flux of atmospheric muon-type neutrinos measured at SuperKamiokande 
PJ, combined with a large deviation of the muon-type to electron- type neu- 
trino flux ratio from theoretical predictions. This "atmospheric neutrino 
puzzle" is best solved by assuming that oscillates into v T and that the v e 
does not oscillate. For a recent analysis of all the atmospheric neutrino data 
see §. 

On the other hand, measurements of the solar neutrino flux |3|, f|, ^ |^. [7j 
have always been plagued by a large suppression of the measured solar v e flux 
with respect to theoretical predictions |J. Again, this "solar neutrino puzzle" 
is best resolved by assuming that v e oscillates into a linear combination of 
the other flavour eigenstates [|9], [K| (for a more conservative analysis of the 
event rates and the inclusion of the "dark side" of the parameter space, see 



II]). The most recent analysis of the solar neutrino data which includes the 



mixing of three active neutrino species can be found in |12 



Neutrino oscillations were first hypothesised by Bruno Pontecorvo in the 



1950's fI3|. The hypothesis of three flavour mixing was first raised by Maki, 



Nakagawa and Sakata |14| . In light of the solar neutrino puzzle, Wolfenstein 



and Mikheyev and Smirnov ||TB[ realized that neutrino-matter interac- 
tions could affect in very radical ways the survival probability of electron- 
type neutrinos which are produced in the solar core and detected at the Earth 
(MSW effect). 

Since then, significant effort has been devoted to understanding the os- 
cillation probabilities of electron-type neutrinos produced in the Sun. For 
example, in [0 the survival probability of solar electron-type neutrinos was 
discussed in the context of three-neutrino mixing including matter effects, 
and solutions to the solar neutrino puzzle in this context were studied (for 



example, in [O, IT 



In this paper, the understanding of solar neutrino oscillations is extend to 
the case of other active neutrino species [y^ z/ T , and antineutrinos) produced 
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in the solar core. Even though only electron-type neutrinos are produced 
by the nuclear reactions which take place in the Sun's innards, it is well 
know that, in a number of dark matter models, dark matter particles can 
be trapped gravitationally inside the Sun, and that the annihilation of these 
should yield a flux of high energy neutrinos (E v ^ 1 GeV) of all species which 
may be detectable at the Earth [[HJ. Indeed, this is one of the goals of very 
large "neutrino telescopes," such as AMANDA |(| or BAIKAL |2l|. It is 
important to understand how neutrino oscillations will affect the expected 
event rates at these experiment s.[] 

The oscillation probability of all neutrino species has, of course, been 
studied in different contexts, such as in the case of neutrinos produced in the 
core of supernovae KBJ or in the case of neutrinos propagating in constant 



electron number densities \\24\\ . The latter case has been receiving a consider- 



able amount of attention from neutrino factory studies [EH] . The case at hand 



(GeV solar neutrinos) differs significantly from these mentioned above, in at 
least a few of the following: source-detector distance, electron number density 
average value and position dependency, energy average value and spectrum. 
Neutrino factory studies, for example, are interested in 0(1000) km base- 
lines, O(10) GeV electron-type and muon-type neutrinos produced via muon 
decay propagating in roughly constant, Earth-like (matter densities around 
3 g/cm 3 ) electron number densities. 

The paper is organised as follows. In Sec. 2, the well known case of two- 
flavour oscillations is reviewed in some detail, while special attention will be 
paid to neutrinos produced inside the Sun. In Sec. 3 the same discussion is 
extended to the less familiar case of three-flavour oscillations. Again, special 
attention is paid to neutrinos produced in the Sun's core. In Sec. 4 the results 
presented in Sec. 3 will be analysed numerically, and the three-neutrino multi- 
dimensional parameter space will be explored. Sec. 5 contains a summary of 
the results and the conclusions. 

It is important to comment at this point that one of the big challenges of 
studying three-flavour oscillations is the multi-dimensional parameter space, 
composed of three mixing angles, two mass-squared differences, and one 
complex phase, plus the neutrino energy. For this reason, the discussions 
presented here will take advantage of the current experimental situation to 
constrain the parameter space, and of the possibility of producing neutrinos 
of all species via dark matter annihilations to constrain the neutrino energies 



"Some effects have already been studied, in the two- neutrino case, in [E2[ 
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to the range from a few to tens of GeV. 



2 Two-Flavour Oscillations 

In this section, the well studied case of two-flavour oscillations will be re- 
viewed [2f|. This is done in order to present the formalism which will be 



later extended to the case of three-flavour oscillations and describe general 
properties of neutrino oscillations and of neutrinos produced in the Sun's 
core. 



2.1 Generalities 

Neutrino oscillations take place because, similar to what happens in the 
quark sector, neutrino weak eigenstates are different from neutrino mass 
eigenstates. The two sets are related by a unitary matrix, which is, in the 
case of two-flavour mixing, parametrised by one mixing angle i9.f\ 

f Ve \_ f U e i U e2 \ / V X \ f COS SU1 ■& \ / V X \ . . 

UJ \U xl U x2 )\v 2 ) v-sin^ costfJW' { -> 

where v\ and v 2 are neutrino mass eigenstates with masses m\ and m 2 , re- 
spectively, and v x is the flavour eigenstate orthogonal to u e . All physically 
distinguishable situations can be obtained if0<$<7r/2 and m\ < m\ or 
< •& < 7r/4 and no constraint is imposed on the masses-squared. 

In the case of oscillations in vacuum, it is trivial to compute the proba- 
bility that a neutrino produced in a flavour state a is detected as a neutrino 
of flavour (3, assuming that the neutrinos are ultrarelativistic and propagate 
with energy E v : 

P*/3 = \U^\ 2 \U al \ 2 + \Up 2 \ 2 \U a2 \ 2 + 2Re (u^U^U^e^^j . (2.2) 

Here Am 2 = m\ — m\ is the mass-squared difference between the two mass 
eigenstates and x is the distance from the detector to the source. It is trivial 
to note that P a/ g = Pp a since all U a i are real and the theory is T-conserving. 
Furthermore, note that d is indistinguishable from tt/2 — 19 (or, equivalently, 

tlf the neutrinos are Majorana particles, there is also a Majorana phase, which will be 
ignored throughout since it plays no role in the physics of neutrino oscillations. 
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the sign of Am 2 is not physical), and all physically distinguishable situations 
are obtained by allowing < $ < tt/A and choosing a fixed sign for Am 2 . 

In the case of nontrivial neutrino-medium interactions, the computation 
of can be rather involved. Assuming that the neutrino-medium inter- 
actions can be expressed in terms of an effective potential for the neutrino 
propagation, one has to solve 



d. /ux(t) 
dt W*) 



—i 



Ex \ / V u (t) V 12 (t) 
E 2 ) + \V l2 (ty V 22 {t) 



U(t) 



(2.3) 



with the appropriate boundary conditions (either a v e or a u x as the ini- 
tial state, for example). In the ultrarelativistic limit one may approximate 
E 2 — E x ~ Am 2 /2E U , d/dt =~ d/dx, and Vij(t) ~ Vij(x). A very cru- 
cial assumption is that there is no kind of neutrino absorption due to the 
neutrino-medium interaction, i.e., the 2x2 Hamiltonian for the neutrino 
system is Hermitian. 

It is interesting to argue what can be said about P a/ 3 in very general 
terms. First, the conservation of probability requires that 

P ee + P ex = 1, (2.4) 

Pxe 4" Pxx = 1- (2-5) 

Second, given that the Hamiltonian evolution is unitary, 

P ee + P xe = 1. (2.6) 

It is easy to show that the extra constraint P ex + P xx = 1 is redundant. 
Eq. (|2.6|) can be understood by the following "intuitive" argument: if the 
same amount of v e and v x is produced, independent of what happens to v e 
and v x during flight, the number of v e and v x detected in the end has to 
be the same. In light of the constraints above, one can show that there is 
only one independent P a p, which is normally chosen to be P ee . The others 
are given by P ex = P xe = 1 — P ee and P xx = P ee . Note that the equality 
P ex = P xe is not a consequence of T-invariance, but a consequence of the 
unitarity of the Hamiltonian evolution and particular only to the two-flavour 
oscillation will be shown later. 



2.2 Oscillation of Neutrinos Produced in the Sun's Core 

It is well known f|TJ, [T6j that neutrino-Sun interactions affect the oscillation 
probabilities of neutrinos produced in the Sun's core in very nontrivial ways. 
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Indeed, all but one solution to the solar neutrino puzzle rely heavily on 
neutrino-Sun interactions || [LC], P| . The survival probability of electron- 
type solar neutrinos has been computed in many different approximations 
by a number of people over the years, and can be understood in very simple 
terms [p6[] . 

In the presence of electrons, the differential equation satisfied by the two 
neutrino system is, in the flavour basis, 



d 

dx 



—i 



Am 2 



\u e2 



u* e2 u, 2 



,u2 



u, 



H2\ 



+ 



(A{x) 
\ 



(X) 

(x) 



(2.7) 

where terms proportional to the 2x2 identity matrix were neglected, since 
they play no role in the physics of neutrino oscillations. 



A(x) = V2G F N e 



x 



(2.8) 



is the charged current contribution to the v e -e forward scattering ampli- 
tude, Gf is Fermi's constant, and N e (x) is the position dependent electron 
number density. In the case of the Sun (see also p7|), A = A(0) ~ 
6 x 10~ 3 eV 2 /GeV, assuming an average core density of 79 g/cm 3 , and A(x) 
falls roughly exponentially until close to the Sun's edge. It is safe to say that 
significantly far away from the Sun's edge A(x) is zero. 

A particularly simple way of understanding the propagation of electron- 
type neutrinos produced in the Sun's core to the Earth is to start with a 
v e state in the basis of the eigenstates of the Hamiltonian evaluated at the 
production point, |z/ e ) = co s{} m (0) l) + sin {} M (0)\u H ), where \u H ) (|^l)) cor- 
respond to the highest (lowest) instantaneous Hamiltonian eigenstate. The 
matter mixing angle i?m = $m(0) is given by 



cos 2$ 



Am 2 cos 2$ - 2E,A 



M 



(Am 2 ) 2 + A 2 - AE v AAm 2 cos 2$ 



(2.9) 



The evolution of this initial state from the Sun's core is described by an 
arbitrary unitary matrix until the neutrino reaches the Sun's edge. From this 
point on, one can rotate the state to the mass basis and follow the vacuum 
evolution of the state. Therefore, P ee (x), where x is is the distance from the 
Sun's edge to some point far away from the Sun (for example, the Earth), is 



[x 



U 



el 





COS^m 



(2.10) 
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where overall phases in the amplitude have been neglected. The matrix 
parametrised by A, B represents the evolution of the system from the Sun's 
core to vacuum, and also rotates the state into the mass basis.f] Expanding 
Eq. (|2.10 ), and assuming that there is no coherence in the Sun's core between 



Ul and one arrives at the well known expression (these have been first 



derived using a different language in p8| and fl29f) 



I / /\ th, T. 

P ee (x) = Pi cos 2 d + P 2 sin 2 ■& - cos 2$ M JP C (1 - P c ) sin 2$ cos — — + 5 

\ 2E,, 



(2.11) 

where 5 is the phase of AB*, P c = |P| 2 = 1— \A\ 2 is the "level crossing 
probability" , and Pi = 1 — P 2 = | + | (1 — 2P C ) cos 2$ M is interpreted as the 
probability that the neutrino exits the Sun as a V\. 

Eq. (|2.11|) should be valid in all cases of interest, and contains a large 



amount of features. In the case of the solar neutrino puzzle, the neutrino 
energies of interest range between hundreds of keV to ten MeV, and mat- 
ter effects start to play a role for values of Am 2 as high as 10~ 4 eV 2 . In 
the adiabatic limit (P c — > 0) very small values of P ee are attainable when 
cos2-#m — * — 1 and sin 2 $ is small. More generally, in this limit P ee = sin 2 i?. 
This is what happens for all solar neutrino energies in the case of the LOW 
solution^] for solar neutrino energies above a few MeV in the case of the 
LMA solution, and for 400 keV < E v < 1 MeV energies in the case of the 
SMA solution. In the extremely nonadiabatic limit, which is reached when 
Am 2 /2E U <C A, P c — > cos 2 $ and cos2-# M — > —1, the original vacuum oscil- 
lation expression is obtained, up to the "matter phase" 5. This is generically 
what happens in the VAC solution to the solar neutrino puzzle. 

If the electron number density is in fact exponential, one can solve Eq. ( |2.7|) 
exactly [[!(], [31]. For N e (x) = N e (0) e~ x ^ r °, where x = is the centre of the 
Sun, 

-7 sin 2 ■& _ —7 

Pc = —i =^-, (2-12) 



■i-The most general form of a 2 x 2 unitary matrix is I „„ ) f n i C ) > wnere 



-B* A* J \0 e iC 

\A\ 2 + \B\ 2 = 1 and < C < 27r. In the case of neutrino oscillations, however, the physical 
quantities are \ A\ 2 and the phase of AB* , and therefore £ can be ignored. 

§This is in general the case, because one has to consider that neutrinos are produced 
at different points in space and time. 

^See H [l0| [ll[] for the labelling of the regions of the parameter space that solve the 
solar neutrino puzzle 
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where 



7 = 27rr ' 



Am 2 
2E~ 



1.05 



Am 2 



lO- 6 eV z 



1 GeV 



(2.13) 



for ro = -R0/IO.54 = 6.60 x 10 4 km p? |. In the case of the Sun, the expo- 
nential profile approximation has been examined ||32|| , and was shown to be 



very accurate, especially if one allows r to vary as a function of Am 2 /2E U . 

The exact expression for S has also been obtained p§| , and the readers 
are referred to |33| for details concerning physical implications of the matter 
phase. Its effects will not be discussed here any further. 



2.3 The Case of Antineutrinos 

Antineutrinos that are produced in the Sun's core obey a differential equation 
similar to Eq. fl2.7| ), except that the sign of the matter potential changes, i.e. 
A(x) <-> —A(x), and U a i <-> U*^ (this is immaterial since, in the two-flavour 
mixing case, all U a i are real). 

Instead of working out the probability of an electron-type antineutrino 
being detected as an electron-type antineutrino P gg from scratch, there is 
a very simple way of relating it to P ee . One only has to note that, if the 
following transformation is applied to Eq. ( |2.7|) : $ — > n/2 — subtract 
the matrix A(l2 X 2), where I2X2 is the 2x2 identity matrix and relabel 
v e [x) <-> v x [x\ the equation of motion for antineutrinos is obtained.^] There- 
fore, P gg (tf) = P xx (n/2 - ft) = P ee (n/2 - ft) (this was pointed out in 0). 
Remember that, in the case of vacuum oscillations, $ is physically equivalent 
to 7r/2 — so P S e = P e e- In the more general case of nontrivial matter ef- 
fects, this is clearly not the case, since the presence of matter (or antimatter) 
explicitly breaks CP-invariance. 

It is curious to note that, in the case of two-flavour oscillations, there 
is no T-noninvariance, i.e., P a/ 3 = Pp a , while there is potentially large CP 
violation, i.e., P a p 7^ P g ^, even if the Hamiltonian for the system is explicitly 
T-noninvariant and CP-noninvariant, as is the case of the propagation of 
neutrinos produced in the Sun (namely A(t) is a generic function of time and 
A(t) for neutrinos is — A(t) for antineutrinos). 

"if one decides to limit < < 7r/4, a similar result can be obtained if Am 2 — > —Am 2 , 
explicitly P gg (Am 2 ) = P ee (~Am 2 ). 
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3 Three Flavour Oscillations 



Currently, aside from the solar neutrino puzzle, there is an even more convinc- 
ing evidence for neutrino oscillations, namely the suppression of the muon- 
type neutrino flux in atmospheric neutrino experiments [|l|. This atmospheric 
neutrino puzzle is best solved by <-> v r oscillations with a large mixing 
angle 0. Furthermore, the values of Am 2 required to solve the atmospheric 
neutrino puzzle are at least one order of magnitude higher than the values 
required to solve the solar neutrino puzzle. For this reason, in order to solve 
both neutrino puzzles in terms of neutrino oscillations, three neutrino families 
are required. 

In this section, the oscillations of three neutrino flavours will considered. 
In order to simplify the discussion, I will concentrate on neutrinos with ener- 
gies ranging from a few to tens of GeV, which is the energy range expected 
for neutrinos produced by the annihilation of dark matter particles which 
are possibly trapped inside the Sun. Furthermore, a number of experimen- 
tally inspired constraints on the neutrino oscillation parameter space will be 
imposed, as will become clear later. 



3.1 Generalities 

Similar to the two-flavour case, the "mapping" between the flavour eigen- 
states, u e , and v T and the mass eigenstates z/j, % — 1,2, 3 with masses m ; 
can be performed with a general 3x3 unitary matrix, which is parametrised 
by three mixing angles (9, to, and £) and a complex phase 0. In short hand 
notation v a = Ua.i v i where a — e, fi, r and i — 1, 2, 3. The MNS mixing ma- 
trix will be written, similar to the standard CKM quark mixing matrix 
35 , as 



U e i U e2 U e3 \ I cuj c£ suj c£ s£e' 

U,t\ U^ 2 U^\ = \ -suj c9 - cuj s9 s£e _i * cuj c9 - slo s9 sfe - ** s9 c£ 
U Tl U r2 U t3 J \ su> s9 — cuj c9 ste-** -cuj s9 - suj c9 e -i * c9 c£ , 

(3.1) 

where cQ = cos^ and s( = sin( for ( = u,9,£. If the neutrinos are 
Majorana particles, two extra phases should be added to the MNS matrix, 
but, since they play no role in the physics of neutrino oscillations, they can 
be safely ignored. All physically distinguishable situations can be obtained if 
one allows < <fi < tt, all angles to vary between and 7r/2 and no restriction 
is imposed on the sign of the mass-squared differences, Am-. = m 2 — m 2 . 
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Note that there are only two independent mass-squared differences, which 
are chosen here to be Am^ and Am^. 

All experimental evidence from solar, atmospheric, and reactor neutrino 
experiments [fH, j3|, fl[ p|, [7|, p6| can be satisfied,^ somewhat conservatively, by 
assuming |26|: 10~ 4 eV 2 < \Am 2 31 \ ~ |Am 2 2 | < 10~ 2 eV 2 , 0.3 < sin 2 9 < 0.7, 
10~ n eV 2 < lAmfJ ^ 10~ 4 eV 2 , sin 2 £ < 0.1, while is mostly uncon- 
strained. There is presently no information on 0. In determining these 
bounds, it was explicitly assumed that only three active neutrinos exist. 

A few comments about the constraints imposed above are in order. First, 
one may complain that u is more constrained than mentioned above by 
the solar neutrino data. The situation is far from definitive, however. As 
pointed out recently in [jD]] if the uncertainty on the 8 B neutrino flux is 
inflated or if some of the experimental data is not considered (especially the 
Homestake data 0) in the fit, a much larger range of Am^ and u is allowed. 



Furthermore, if three-flavour mixing is considered ||12| , different regions in the 
parameter space Am^-sin 2 uj are allowed for different values of sin 2 £, even 
if sin 2 £ is constrained to be small. 



Second, the limit from the Chooz and Palo Verde reactor experiments |36 
do not constrain sin 2 £ for jAm^l 10~ 3 eV 2 . Furthermore, their constraints 
are to sin 2 2£, so values of sin 2 £ close to one should also be allowed. However, 
the constraints from the atmospheric neutrino data require cos 2 £ to be close 
to one. This is easy to understand. Assuming that L^f is much larger than 
the Earth's diameter and that Am 2 x = Am| 2 , 



P^ m = 1-4 cos 2 £ sin 2 9(1 - cos 2 £ sin 2 9) sin 2 ( =fg±^ J , (3.2) 



Am 2 ] a; 

according to upcoming Eq. ( |3.3| ) . Almost maximal mixing implies that 
cos 2 £ sin 2 # ~ 1/2. With the further constraint from P^ m , namely sin 2 2£ ~ 
0, one concludes that cos 2 £ ~ 1 and sin 2 9 ~ 1/2. 

In the case of oscillations in vacuum, it is straight forward to compute 
the oscillation probabilities P a p of detecting a flavour (3 given that a flavour 
a was produced. 

Am? . x 

p a p = E u ai u aj u pi u;/^- (3.3) 

hi 



* There is evidence for neutrino oscillations coming from the LSND experiment |]37| . 
Such evidence has not yet been confirmed by another experiment, and will not be con- 
sidered in this paper. If, however, it is indeed confirmed, it is quite likely that a fourth, 
sterile, neutrino will have to be introduced into the picture. 
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The three different oscillation lengths, L l J sc , are numerically given by 



L% c = ^ = 2-47 x 10 8 km ( ^ ) ( ) , (3.4) 



ij 




which are to be compared to the Earth-Sun distance (1 a.u.= 1.496 x 10 8 km). 
In the energy range of interest, 1 Gev< E u < 100 GeV and given the ex- 
perimental constraints on the parameter space described above, it is easy 
to see that L^ c and L^ c are much smaller than 1 a.u., and that its effects 
should "wash out" due to any realistic neutrino energy spectrum, detector 
energy resolution, or other "physical" effects. Such terms will therefore be 
neglected henceforth. In contrast, L^l c maybe as large as (and maybe even 
much larger than!) the Earth-Sun distance. Note that a nonzero phase <ft 
implies T- violation, i.e., P a p ^ Pp a , unless L^ c ^> 1 a.u.. This will be 
discussed in more detail later. 

In the presence of neutrino-medium interactions, the situation is, in gen- 
eral, more complicated (indeed, much more!). Similar to the two-neutrino 
case, it is important to discuss what is known about the oscillation probabil- 
ities. From the conservation of probability one has 

Pee Pefi Per 1; 
p _|_ P J_ p — 1 

and, similar to the two-neutrino case, unitarity of the Hamiltonian evolution 
implies 

p _L p I p — i 

Pev + Pw+Prn = ( 3 - 6 ) 

A third equation of this kind, P er + P^ T + P TT = 1, is redundant. As before, 
Eqs. ( |3.6j ) can be understood by arguing that, if equal numbers of all neutrino 
species are produced, the number of z/g's to be detected should be the same, 
regardless of /3, simply because the neutrino propagation is governed by a 
unitary operator. 

One may therefore express all P a p in terms of only four quantities. Here, 
these are chosen to be P ee , P eA1 , P^, and P TT . The others are given by 

P = 1 - P - P 

1 er 1 1 ee 1 ejii 
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p = 1 + p -P -P -P 

1 fie x ~ 1 tt 1 ee - 1 fif- - 1 e^; 

-f/ir -fee Pen P-r; (^■'O 

P = P + P — P 

P = 1 - P - P 



Note that, in general, P a p 7^ P/ 



/3a- 



3.2 Oscillation of Neutrinos Produced in the Sun's Core 



The propagation of neutrinos in the Sun's core can, similar to the two- 
neutrino case, be described by the differential equation 



dx 



(3.8) 



where 5 V $ is the Kronecker delta symbol. Terms proportional to the iden- 
tity 8 a p are neglected because they play no role in the physics of neutrino 
oscillations. The matter induced potential A(x) is given by Eq.( |2.8|) . 

As in the two-neutrino case, it is useful to first discuss the initial states v a 
in the Sun's core, and to express them in the basis of instantaneous Hamilto- 
nian eigenstates, which will be referred to as \vh), \ v m)i an d (H = high, 
M= medium, and L— low). Therefore 



\v a ) = H a \v H ) + M a \v M ) + L a \u L ), 



(3.9) 



where (v a \v a i) = 5 aa >. As before (see Eq. |2.10| ), the probability of detecting 
this initial state as a /3-type neutrino far away from the Sun (e.g., at the 
Earth) is given by 



a/3 



(I 



Upi Up 2 U£ 3 







e 1 2E » 
\0 












(^3x3) 




/ 





(3.10) 

where 14 X 3 is an arbitrary 3x3 unitary matrix which takes care of propa- 
gating the initial state until the edge of the Sun and rotating the state into 
the mass basis. 

In order to proceed, it is useful take advantage of the constraints on 
the neutrino parameter space and the energy range of interest. Note that 
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A > ^™ 31 3> (remember that the energy range of interest is 1 GeV< 

E v < 100 GeV and that A ~ 6 x 10 -3 eV 2 / GeV). It has been shown explicitly 
TTfl , assuming the neutrino mass-squared hierarchy to be m| > m 2 , > mf, [] 



that, if the mass-squared differences are very hierarchical (lAmijjJ 3> jAm 2 ,^), 
the three-level system "decouples" into two two- level systems, i.e., one can 
first deal with matter effects in the U H—M" system and then with the matter 
effects in the "M — L" system. One way of understanding why this is the 
case is to realize that the "resonance point" corresponding to the Am^ is 
very far away from the resonance point corresponding to Amf^ . With this in 
mind, it is fair to approximate (this is similar to what is done, for example, 
in 



38 





f 


A L 


B L 




(i 








^3x3 = 




-B L 




:) 





A H 


B H 




{ 










v0 


_ B H* 


A H * 


where \B H \ 2 = 1 - \A 


H 


2 — 


r c i \ D \ 


1 - 


\A L \ 2 = 


P C L . 



(3.11) 



The superscripts 

H, L correspond to the "high" and the "low" resonances, respectively. 

It also possible to obtain an approximate expression for the initial states 
in the Sun's core. Following the result outline above, this state should be 
described by two matter angles, £m and u M , corresponding to each of the 
two- level systems. Both should be given by Eq. (|2.9|), where, in the case of 



cos2£ M , $ is to be replaced by £ and Am 2 by Am^, while in the case of 
cos2o;m, $ is to be replaced by u>, Am 2 by Am^ and A is to be replaced by 
A cos£ ]38], [L^]. Furthermore, because A cos£ ^> ^^ 21 ^ , cos 2u>m can be safely 
replaced by -1 (remember that cos 2 £ ^ 0.9). Within these approximations, 
in the Sun's core, 

\u e ) = sin £ M | p H ) + cos £ M \u M ), 

\u^) = sin 9 cos $,mWh) — siia 9 sm^ M \u M ) — cos OI^l), (3.12) 
\v T ) = cos9cos£m\vh) —cos 9 sin ^ m \vm) + sin- 9\i>l)- 

The accuracy of this approximation has been tested numerically in the range 
of parameters of interest, and the difference between the "exact" result and 
the approximate result presented in Eq. (|3.12|) is negligible. 

Keeping all this in mind, it is straight forward to compute all oscillation 
probabilities, starting from Eq. ( |3.10| ). From here on, = (no T- violating 



phase in the mixing matrix, such that all U a i are real) will be assumed, in 
n will work under this assumption for the time being. 
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order to simplify expressions and render the results cleaner. In the end of 
the day one obtains 

( /\ 777, 07 \ 

or (3.13) 
Pa/3 = (a a U pi + b a U p2 f + c 2 (Uf^Y - 4a a b a (U pi U p2 ) sin 2 {^^- + 6 L ^j , 

where 5 L is the matter phase, induced in the low resonance, and 



P2 H P C L , 



be = JPfil-P}), 



I pH 
^3 ) 



a„ = - J ( 1 - P C L ) cos 9 - JP^PL sin 9 



^ = yjpr. cos e- tJP?(1 -P c l ) sin^, (3.14) 



-I-' 



Pn H sin 2 9, 



a T = J {I -P C L ) sin 9- JP^P^ cos 9, 



b T = -JP^sm9- J P^ (1 - P^) cos 9, 



and P«f = 1 - P 3 H = (\A H \ 2 cos 2 ^ M + \B H \ 2 sin 2 f M ), which can also be 
written as P^ = \ + \ (l — 2P C H ) cos 2£m- This is to be compared with the 
expression for P\ obtained in the two-flavour case. Note that a 2 +6 2 +c 2 = 1. 
The effect of 5 L will not be discussed here and from here on 5 L will be set 
to zero. For details about the significance of S L for solar neutrinos in the 
two-flavour case, readers are referred to |3^| . 

Many comments are in order. First, in the nonadiabatic limit which 
can be obtained for very large energies, P^ — ► cos 2 £, P 6 L — > cos 2 to and 
cos2£m — ► — 1. It is trivial to check that in this limit a a — > 6 a — > {7 a 2, 
c a — > U a3 , and the vacuum oscillation result is reproduced, up to the matter 
induced phase 5 L . 

Second, P ee can be written as 

P ee = P 2 H cos 2 ^Pln + P S H sm 2 e, (3-15) 
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where P 2 ^ is the two-neutrino result obtained in the previous section (see 
Eq. ( |2.11| )) in the limit cos2$ M — > — 1. It is easy to check that Eq. (|2.11|) 
would be exactly reproduced (with $m replaced by u>m, of course) if the 
cos2ujm = —1 approximation were dropped. 

For solar neutrino energies (100 keV< E v < 10 MeV), £ M -> C, -> 
and therefore P^ (P-i 1 ) cos 2 £ (sin 2 £), reproducing correctly the result 
of the survival probability of electron-type solar neutrinos in a three-flavour 
oscillation scenario (see p6[ and references therein). In this scenario there is 



no U H — L resonance inside the Sun, because 2E 31 ^> A for solar neutrino 
energies. 

On the other hand, in the case P^ — > and cos 2^/ —* — 1, P$ 1 
and electron-type neutrinos exit the Sun as a pure u 3 mass eigenstate, and 
do not undergo vacuum oscillations even if Am^ is very small. In contrast, 

and v T always undergo vacuum oscillations if Am^ is small enough. The 
reason for this is simple. The generic feature of matter effects is to "push" 
v e into the heavy mass eigenstate, while and u T are "pushed" into the 
light mass eigenstates. This situation is changed by nonadiabatic effects, as 
argued above. 

Finally, it is important to note that all equations obtained are also valid 
in the case of inverted hierarchies (m 2 < m 2 2 or m 2 < m 2 ). This has 
been discussed in detail in the two- neutrino oscillation case [^9|, and is also 
applicable here. It is worthwhile to point out that, in the approximation 
Awijj ~ Am| 2 the transformation Am^i — > — Am|i can be reproduced by 
transforming u — ► tt/2—uj, 9 — > n — 6 and redefining the sign of v T . Therefore, 
one is in principle allowed to fix the sign of A long as 6 is allowed to 

vary between and ir. 

In the case of inverted hierarchies (especially when Am 2 x < 0) one expects 
to see no "level crossing" (indeed, matter effects tend to increase the distance 
between the "energy" levels in this case), but matter effects are still present, 
because the initial state in the Sun's core can be nontrivial ( i.e., $m 7^ 
Note that v e is still "pushed" towards Vh, even in the case of inverted hi- 
erarchies, and the expressions for the matter mixing angles Eq. ( |2.9|) , and 
the initial states inside the Sun Eq. ( |3.12|) are still valid. The consequence 
of no "level crossing" is that the adiabatic limit does not connect, for exam- 
ple, Uh —> v 3 but Vh — > v 2 (or V\, depending on the sign of Am^). This 
information is in fact contained in the equations above. The crucial feature 
is that, for example, when Am 2 x < 0, P^ 1 — > 1 in the "adiabatic limit," 
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and the matrix V3X3 correctly "connects" uh — > ^2 (or V\)\ Another curious 
feature is that, in the limit \Am^ x \/2E v ^> A, cos2£ A / — > — cos 2^, P^ — >■ 1 
and Eq. ( |3.15| ) correctly reproduces the survival probability of electron-type 
solar neutrinos in the three-flavour oscillation case. Note that on this case 
the sign of Am^ does not play any role, as expected. On the other hand, it 
is still true that P^ H,L ' — > cos 2 (£, a;) in the extreme nonadiabatic limit, and 
vacuum oscillation results are reproduced, as expected. Again, in this limit, 
one is not sensitive to the sign of Am^, as expected. 



3.3 The Case of Antineutrinos 

As in the two-neutrino case, the difference between neutrinos and antineutri- 
nos is that the equivalent of Eq. ( p.8|) for antineutrinos can be obtained by 
changing A(x) — > —A(x) and U ai <-» U* { . Unlike the two-flavour case, how- 
ever, there is no set of variable transformations that allows one to exactly 
relate the differential equation for the neutrino and antineutrino systems. 
One should, however, note that if the signs of both Am 2 are changed and 
U a i <-> U*i, the neutrino equation turns into the antineutrino equation, up 
to an overall sign. This means, for example, that the instantaneous eigenval- 
ues of the antineutrino Hamiltonian can be read from the eigenvalues of the 
neutrino Hamiltonian with Amfj <-> -Am|, U a i *-> plus an overall sign. 

When it comes to computing P a g this global sign difference is not relevant, 
and therefore P a/ g(Ara|, U ai ) = P a(S (-Am|, C/*J. 



4 Results and Discussions 

This section contains the compilation and discussion of a number of results 
concerning the oscillation of GeV neutrinos of all species produced in the 
Sun's core. The goal here is to explore the multidimensional parameter space 
spanned by Am^, Am^, sin 2 u;, sin 2 9, and sin 2 £ (and E v ). 

It will be assumed throughout that the electron number density profile of 



the Sun is exponential, so that Eq. ( 2.12 ) can be used. As mentioned before, 



the numerical accuracy of this approximation is quite good, and certainly 
good enough for the purposes of this paper. Therefore, both P^ and P C L 
which appear in Eq. (|3.13|) will be given by Eqs. (|2.12| , |2.13|) , with $ — > £, 



Am 2 — > Am^ in the former, and d — > u, Am 2 — > Am 2 ^ in the latter. 

When computing P a/3 , an averaging over "seasons" is performed, which 
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"washes out" the effect of very small oscillation wavelengths. Furthermore, 
integration over neutrino energy distributions is performed. Finally, all P a p 
to be computed should be understood as the value of P a p in the Earth's 
surface, i.e., Earth matter effects are not included. This is done in order to 
make the Sun matter effects in the evaluation of P a p more clear. It should be 
stressed that Earth matter effects may play a significant role for particular 
regions of the parameter space, but the discussion of such effects will be left 
for another opportunity. 

Because the parameter space to be explored is multidimensional, it is 
necessary to make two-dimensional projections of it, such that "illustrative" 
points are required. The following points in the parameter space are chosen, 
all inspired by the current experimental situation: 



• 


ATM: Am 2 31 


= 3 x 10" 


" 3 eV 2 , 


sin 2 9 ■ 


= 0.5, and sin 2 £ = 0.01, 


• 


LMA: Am\ x 


= 2 x 10- 


" 5 eV 2 , 


sin 2 a; 


= 0.2, 


• 


SMA: Am 2 21 


= 6 x 10" 


" 6 eV 2 , 


sin 2 uj 


= 0.001, 


• 


LOW: Am 2 21 


= 1 x 10 


" 7 eV 2 , 


sin 2 a; 


= 0.4, 


• 


VAC: Am§! 


= 1 x 10- 


10 eV 2 , 


sin 2 uj 


= 0.55. 



ATM corresponds to the best fit point of the solution to the atmospheric 
neutrino puzzle ||, and a value of sin 2 £ = 0.01 which is consistent with all 
the experimental bounds. Note that some "subset" of ATM will always be 
assumed (for example, sin 2 9 is fixed while exploring the (Am^ x sin 2 £)- 
plane). For each analysis, it will be clear what are the "variables" and what 
quantities are held fixed at their "preferred point" values. All other points 
refer to sample points in the regions which best solve the solar neutrino puzzle 
H [10], [TT[] , and the notation should be obvious. Initially a flat neutrino 
energy distribution with E™ m = 1 GeV and E™ ax = 5 GeV is considered 
(for concreteness) , and the case of higher average energies is briefly discussed 
later. 

4.1 The Case of Vacuum Oscillations 

If neutrinos were produced and propagated exclusively in vacuum, the oscil- 
lation probabilities would be given by Eq. ( |3.3| ). This would be the case of 
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neutrinos produced in the Sun's core if either the electron number density 
were much smaller than its real value or if very low energy neutrinos were 
being considered. Nonetheless, it is still useful to digress some on the "would 
be" vacuum oscillation probabilities in order to understand better the matter 
effects. 

In the case of pure vacuum oscillations, it is trivial to check that P a p = 
Pp a (remember that the MNS matrix phase cf) has been set to zero), and 
therefore all P a $ can be parametrised by three quantities, namely P aa , a = 
e, /i, r. It is easy to show that 

P aP = & Pep. = + Prr ~ ~ Pee)- (4.1) 



From Eq. (J^ 

Paa = ^l + ^ 4 2 + ^ 4 3 + 2^ 2 COs(^i) 

or (4.2) 
Paa = (l-^ 3 ) 2 + ^3-4^ 2 sin 2 (^). 

Note that there is no dependency on Am^. Particularly simple limits can 
be reached when L 2 g C is either very small or very large compared with the 
Earth-Sun distance. In both limits P aa is independent of Am^ and, in the 
latter case, P aa depends only on [/ 2 3 . Fig. [I] depicts constant P a p contours 
in the (Am^ x sin 2 u;)-plane, at ATM. Remember that, here, P efl is not an 
independent quantity but is a linear combination of all P aa . Note that P ee 
is symmetric for uj — > n/2 — oj, and that P w <-» P TT when oj — > tt/2 — uj. 
The latter property is a consequence of 6 = ir/4. Also, in the case of P ee , 
the LH C — > oo coincides with the uj — > 0,7r/2 limit for any LH C (this is 
because either U e \ or U e 2 go to zero). This is not true of P MM or P TT unless 
sin 2 ^ = 0. Another important consequence of L 2 g C ^> 1 a.u. is that T- 
violating effects are absent, even if is nonzero. This can be seen by looking 
at the second expression in Eq. ( |4.2|) , which is a function only of |t/ a 3| 2 in 
the limit L 2 g C — ► oo. 

Finally, one should note that oscillatory effects are maximal for Am^ ~ 
2 x 1CT 8 eV 2 . In this region cos (■^^ J — ) — — 1, and the largest suppression to 
all P aa is obtained when U^U^ is maximum. For example, P ee is smallest 
when uj = 7r/4, since U^U^ oc sin 2 2u. There are no "localised" maxima for 
P aa because ^ 2 !^ 2 2 is positive definite. 
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Figure 1: Constant P a p contours in the (Am^ x sin 2 cj)-plane, at ATM (see 
text), in the case of pure vacuum oscillations. 

4.2 "Normal" Neutrino Hierarchy 

When matter effects are "turned on," the situation can be dramatically dif- 
ferent. This is especially true in the case of normal neutrino mass hierarchies 
(m\ < m\ < m 2 ), which will be discussed first. 

The first effect one should observe is that, even though LH C <C 1 a.u., 
Pap depend rather nontrivially on Am^. This dependency comes from the 
terms P 3 H and P^ = 1 — Pf in Eq. (|3.13|) . Remember that P^ is interpreted 
as the probability that a u e produced in the Sun's core exits the Sun as a 
u 3 mass eigenstate. When matter effects are negligible (such as in the limit 
of small neutrino energies) P% — > sin 2 £, its "vacuum limit." Fig. |^ depicts 
constant P% contours in the (Am^/ E u x sin 2 £)-plane. 

Note that, for Am\ x jE v ~ 1CT 2 eV 2 /GeV, P% — > 1, even for small values 
of sin 2 £. In this region, z^'s produced in the Sun's core exit the Sun as pure 
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z/ 3 's. Therefore, P ea ~ £7 2 3 . Because of unitarity in the propagation, v^s and 
Vt's exit the Sun as linear combinations of the light mass eigenstates, and 
may not only undergo vacuum oscillations but are also susceptible to further 
matter effects (dictated by the "M — L" system, as described in Sec. 3). For 
future reference, at ATM, P% ~ 0.87 when averaged over the energy range 
mentioned in the beginning of this section. 

As Aml 1 /E u decreases (as is the case for higher energy neutrinos) the 
nonadiabaticity of the U H — M" system starts to become relevant, and P 3 H — > 
sin 2 £, as argued in Sec. 3.2. A hint of this behaviour can already be seen in 
Fig. [|, for small values of Am^/E,,. 

The information due to the 11 M — V matter effect is encoded in P C L , 
present in Eq. ( |3.13|) . Fig. || depicts contours of constant 1 — P C L in the 
(Am,2i/E u x sin 2 u;)-plane. One should note that 1 — P C L reaches its extreme 
nonadiabatic limit, sin 2 u, when Am^/Pj, < 10~ 7 eV 2 /GeV. For Am| 1 /-E' I/ ^ 
1CT 7 eV 2 /GeV, matter effects increase the value of 1 — P C L . 

One can use the intuition from the two-flavour solution to the solar neu- 
trino puzzle to better appreciate the results presented here. In the case of 
the solutions to the solar neutrino puzzle, the energies of interest range from 
100 keV to 10 MeV, and large matter effects happen around Am 2 ~ 10 -5 eV 2 . 
Furthermore, at Am 2 ~ 10~ 10 eV 2 one encounters the "just-so" solution, 
which is characterised by very long wave-length vacuum oscillations. Rescal- 
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Figure 3: Contours of constant 1 — P C L (see text) in the (Am^/E^ x sin 2 em- 
plane. 

ing to O(GeV) energies, the equivalent of the "just-so" solution happens for 
Amji ~ (10~ 8 — 10 -7 ) eV 2 , while large matter effects would be present at 
Am 2 ~ (10~ 3 — 10~ 2 ) eV 2 . Indeed, one observes large matter effects for 
Am§! ~ (10~ 3 - 10~ 2 ) eV 2 . Am^ ~ (10~ 5 - KT 6 ) eV 2 corresponds to the 
region between the LOW and VAC solutions, where matter effects distort 
P a/ 3 from its pure vacuum value, but no dramatic suppression or enhance- 
ment takes place. Incidently, this behaviour has physical consequences in the 
solution to the solar neutrino problem, as was first pointed out in ||40| . 

Figs. [4] and |5] depict contours of constant P aa and P eAt in the (Am^ x 
sin 2 £)-plane. As expected, in the region where P% ~ 1, P ee and P eil do 
not depend on Am^ or sin 2 u;, namely P ee ~ sin 2 £ and P e}l ~ 0.5 cos 2 £. 
Remember that the results depicted in Figs. |] and | (and all other plots 
from here on) are for an energy band from 1 to 5 GeV. On the other hand, 
P^n and P TT do depend on the point (LMA, SMA, etc), even for P^ ~ 1, as 
foreseen. This dependence will be discussed in what follows. 

In the limit P^ = 1, sin 2 9 = 0.5 

<= 4= 0, (4.3) 
a\= h\= 0.5(1 + 2^(1-^)), (4.4) 
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b]= a 2 T = 0.5(1-2^(1-^)), (4.5) 
aA= -ar&r= 0.5(1 -2P C L ), (4.6) 

and 



± (l-2P c ^)?7 (M)T)1 [/ (/i , r)2 cos 



2P„ 



At both LMA and SMA, the oscillatory term averages out to zero, while 
at VAC cos (-^f 1 -) = 1- It is only at LOW that the oscillatory term is 
nontrivial, as was mentioned in the analogy between the situation at hand 
and the solutions to the solar neutrino puzzle. 

Furthermore, at SMA, 1 — P C L is tiny (see Fig. |), so it is fair to ap- 
proximate P MM ~ P TT ~ 0.5(1 — 0.5 x 0.99) ~ 0.25, in agreement with 
Fig. U(bottom). At LMA, it is fair to approximate sin 2 £ = 0. In this limit, 
U 2 X — U 2 2 ~ — 0.5cos2u; + sin2co>sin£, while U%\ — U% 2 — — 0.5cos2u; — 
sin2c<jsin£. Therefore, because cos2u; = 0.6 > 0, P w is significantly less 



than P TT , since yP c L (l — P C L ) is nonnegligible. Roughly, P m ~ 0.15 and 
Pufi — 0.4, using the approximations above. Again, there is agreement with 
Fig.|(top). 

In order to understand the behaviour at LOW and VAC, one should take 
advantage of the fact that 1 — P C L — > sin 2 uj. In this case, it proves more 
advantageous to use the second form of Eq. ( |3.13|) in order to express all P a p 



I eel 



P ee = P 2 H cos 2 «£ + P 3 H sin 2 £- (Osc) 
P e , = (P 2 H sin 2 e + P 3 ^cos 2 e)sin 2 0-(Osc) eM , (4.8) 

P w = (cos 2 fl + a/p^ sin£sin 2 ^ + P 2 H cos 2 f sin 4 9 - (Osc) 



rr ' 



P TT = (sin 2 6 + ^JPi 1 sin £ cos 2 6 j + cos 2 £ cos 4 6 - [Osc] 
where 

-fif-j ( 4 - 9 ) 

are the oscillatory terms. When L 2 l c ^> 1 a.u., the oscillatory terms are zero, 
and P aj 3 are particularly simple. Note that on this limit many simplifications 
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happen: P a/3 is independent of uj and Am| x , and P w = P TT if sin 2 = cos 2 9, 
as can be observed in Fig. ^(bottom). A very important fact is that, when the 
oscillatory terms are neglected, 2P efl = 1 + P TT — P^ — P ee , as one may easily 
verify directly As argued before, when this condition is satisfied, P a p = Pp a . 
This is not the case in the presence of nonnegligible oscillation effects or when 
P c cos 2 u;. Both statements are trivial to verify directly. For example, 

Aa e b e U^iU^ 2 = Pf sin2u; 
while 



sin 2uj [sin 2 £ sin 2 9 — cos 2 9 J — sin £ sin 29 cos 2uj 

(4.10) 



4a /U 6 At C/ e iC/ e 2 = cos 2 £ sin 2a; 



sin 2uj (P 3 H sin 2 9 - cos 2 0) - yj P" sin 26 cos 2uo 

(4.11) 

SO 0SC 6M ^ OSC^. 

Figs. | and [7] depict P a/3 as a function of sin 2 9 at LMA and SMA, and at 
LOW and VAC, respectively. In these figures, all P a p are plotted, in order to 
illustrate that P aP ^ P Pa at LMA, SMA and VAC. Note that at LMA and 
SMA, the difference comes from the fact that P^ 7^ sin 2 £ and P C L 7^ cos 2 co>. 
At LOW, P C L ~ cos 2 u, but P^ 7^ sin 2 £ and nontrivial oscillatory terms 
render P a/3 7^ Pp a . At VAC, even though P^ 7^ sin 2 £, P a/3 = Pp a because 
P C L = cos 2 uj and because "1-2" oscillations don't have "time" to happen. 

From Eqs. (|4.8|) one can roughly understand the dependency of P a p on 
sin 2 9. Obviously P ee does not depend on 9 (by the very form of the MNS ma- 
trix, Eq. ( p. 1|) ), while P efl (P eT ) depends almost exclusively on sin 2 9 (cos 2 9). 
This is guaranteed by the fact that P^ 1 ^> P^ even at LMA and LOW, 
when one expects the interference terms to play a significant role. It is also 
worthwhile to note that, as expected, at VAC and SMA the curves are very 
similar, a behaviour that can be understood from earlier discussions. 

Finally, Fig. [8] depicts constant P a p contours in the (Am^ x sin 2 u;)-plane, 
at ATM. In light of the previous discussions, the shapes and forms can be 
readily understood. First note that the shapes of the constant P ee and P efl 
regions resemble those of the pure vacuum oscillations depicted in Fig. p], with 
two important differences. First, the constant values of the contours are quite 
different. For example, P ee varies from a few percent to less then 15%, while 
in the case of pure vacuum oscillations, P ee varies from 30% to 100%. This 
can be roughly understood numerically by noting that P( ee , e ^) — P^PJ^eii) 
(remember that P^ 1 = 1 — P% ~ 0.13 when averaged over the energy range 
of interest). Second, at high Am 2 1; the regions are distorted. This is due 
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Figure 6: P a p as a function of sin 2 9, at LMA (top) and SMA (bottom) (see 
text). 
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Figure 7: P a p as a function of sin 2 9, at LOW (top) and VAC (bottom) (see 
text). 
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Figure 8: Constant P a p contours in the (Am^ x sin 2 cj)-plane, at ATM (see 
text). 

to nontrivial matter effects in the "M — L" system. Note that the contours 
follow the constant 1 — P C L curves depicted in Fig. ||. 

The P^ and P TT contours are a lot less familiar, and require some more 
discussion. Many features are rather prominent. For example, the plane is 
roughly divided into a sin 2 u; > 0.5 and sin 2 co> < 0.5 structure, and large 
(small) values of P^ (P TT ) are constrained to the sin 2 u > 0.5 half, and vice- 
versa. Also, there is a rough P J1/1 (w) <-> P TT (n/2— to) symmetry in the picture, 
which was present in the pure vacuum case (see Fig. p. This symmetry is 
absent for large values of Am^, similar to what happens in the case of P ee , 
and is due, as mentioned in the previous paragraph, to the fact that P C L is 
significantly different from cos 2 u in this region. 

The other features are also fairly simple to understand, and are all due 
to fact that P% ^> sin 2 £. It is convenient to start the discussion in the limit 
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when L^i 3> 1 a.u. (the very small Am^ region). As was noted before, 
P a p are given by Eq. (|4.8|) where the Osc a p terms vanish. It is therefore 
easy to see that P a p do not depend on uj or Am^ (as mentioned before), 
and furthermore it is trivial to compute the value of P a p given that we are 
at ATM and that P 3 H ~ 0.87. The next curious feature is that there is a 
"band" around sin 2 uj = 1/2 where P MM ~ P MAt (L^ c — > oo). The same is true 
of P TT . This is due to the fact that, around sin 2 a; ~ 1/2, a^b^ and a T b T 
vanish when P% is large. In the limit P% = 1 one can use Eq. Q4.6|) and note 
that indeed both a^b^ and a T b T vanish at P C L = 1/2. However, for values of 
Anij) ^ 10~ 7 eV 2 P C L ~ cos 2 c<j, which explains the band around uj ~ 7r/4. 
Slight distortions are due to the fact that P 3 H ^ 1, and are easily computed 
from the exact expressions. 



if sin 2 £ terms are neglected. The +, — signs are for P w while the — , + signs 
for P TT . It is trivial to verify numerically (if a little tedious) that the P w 
term has a maximum at sin 2 u; ~ 0.1 and a minimum at sin 2 u; ~ 0.8. For 
P TT the maximum (minimum) is at sin 2 a; ~ 0.9(0.2). It is important to 
comment that the minima are negative numbers. On the other hand, from 
Fig. [TJ (as mentioned before) it is easy to see that cos (— ^rg 51 - J is minimum 
for Amjj ~ 2 x 10~ 8 eV 2 (this is where all P aa are maximally suppressed 
in Fig. |l]). Combining both informations, it is simple to understand the 
maxima/minima of P w and P TT at Am^ — 2 x 10~ 8 eV 2 : Minima occur when 
the coefficient is maximum (e.g., at sin 2 u ~ 0.1 for P MM ) while maxima occur 
when the coefficient is minimum (e.g., at sin 2 uj ~ 0.8 for P^). A description 
of what has happened is the following: The matter effects "compress" the 
constant P^ (P TT ) contours from the pure vacuum oscillation case (presented 
in Fig. [I]) to the sin 2 u; < 1/2 (> 1/2) half of the plane, and a new region 
"appears" on the other half. This other region is characterised by negative 
values to the coefficients of the oscillatory terms, which are not attainable in 
the case of pure vacuum oscillations (see Eq. ( |4.2j )). 

At last, the contours in the region where the oscillatory effects average 
out, and P TT are also best understood from Eq. (|4.7|) and the paragraphs 
which follow it, in the limit that cos (— ) — ► 0. It is simple to see, for 



Again in the limit P% = 1 , P C L 
term Eq. (|4.7|) is 



cos 2 uj, the coefficient of the cos 





(4.12) 
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example, that < P TT if cos 2oo > (sin a; > 1/2), while the situation is 
reversed if cos2u; < 0. This is indeed what one observes in Fig. 

4.3 "Inverted" Neutrino Hierarchy 

Here I turn to the case of an "inverted" neutrino hierarchy, namely Am^ < 0. 
Currently, there is no experimental hint as to what the sign of Am^ should 
be, so there is no reason to believe that the "normal" hierarchy is to be 
preferred over the "inverted" hierarchy. Indeed, even from a theoretical/ 
model building point of view, there are no strong reasons for or against a 
particular neutrino mass hierarchy [f|l[ 



The discussion will be restricted to Am^ > for two reasons. First, 
the Am^ < can be approximately read off from the Am^ > case by 
changing u — > tt/2 — lu, as mentioned before. Second, and most important, 
there is some experimental hints as to what is the sign of Am^ |12|, . For 
example, the SMA solution only exists for one sign of Am^, while the LMA 
and LOW solutions prefer one particular sign. Even in the case of VAC there 
is the possibility of obtaining information concerning the sign of Am^ from 
solar neutrino data j4(J . Therefore, the notation introduced in the beginning 
of this section (ATM, SMA, LMA, LOW, VAC) still applies, and one should 
simply remember that here Am^ < 0. 

As advertised, the largest effect of Am^ < is the typical values of Pf . 
From Eq. ( |2.12j ), keeping in mind that here 7 is negative, 

1 _ p -|7|cos 2 £ 

P " = !_ e -M • < 413 > 

where 7 is given by Eq. (|2 . 1 3|) with Am 2 — > Am^. Since I7I ^> 1 (see 
Eq. (ETIJP), Pf = 1 for all values of Am 2 31 and sin £ of interest. Indeed, 



this is true for any value of sin 2 £ as long as Am^ ^> 1CT 6 eV 2 . This is to 
be contrasted to the normal hierarchy case, when there is always some value 
of sin 2 £ (which is a function of Am 2 x ) below which P c deviates significantly 
from its adiabatic limit. 

All of the Am 2 x dependency of P Q/ g is therefore encoded in £m- However, 
in the case Amf^ < it is trivial to show that — 1 < cos 2£m < — cos 2£, where 
the upper bound is reached in the limit IAm^/2^1 3> A (this has been 
mentioned before. The minus sign takes care of the "unorthodox" P c — > 1 
adiabatic limit). Since one is interested in sin 2 £ < 0.1 (— cos2£ < —0.8), the 
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Figure 9: P af3 as a function of sin 2 £, at LMA, SMA, LOW, and VAC (see 
text), for Am| x < 0. 



range for £m is rather limited, and therefore any Am^ effects are bound to 
be very small. Larger Am| x effects are expected for larger sin 2 £. 

In light of this, Fig. |9| depicts P a p and P eA1 as a function of sin 2 £ at the 
various points (LMA, SMA, LOW, VAC), for Am 2 31 = -3 x 10~ 3 eV 2 and 
sin 2 9 = 0.5. 

It is interesting to compare the results presented here with the pure vac- 
uum case. In the limit P^ = 1 



P ee = COS 2 £ 



P C L cos 2 u + (1 - Pe) sin 2 u + JP c {l- P l c ) sin 2u cos ' ' ' 2 1 



2E. 



V 



(4.14) 



30 



while the pure vacuum result in the same region of the parameter space is 



PIT = cos 4 e 



4 4 2 2 / A.TYVo-\L 

cos uj + sin uj + 2 sin a> cos uj cos 

V , 



+ sin 4 £. 



(4.15) 

In the limit P C L = cos 2 a;, the difference 

P ee - PIT = (P£> vac ) cos 2 £ sin 2 £ - sin 4 f , (4. 16) 

where P^ ,vac is the electron neutrino survival probability in the two-flavour 
case with Am 2 = Am^ and vacuum mixing angle uj. This difference vanishes 
at sin £ = 0, and sin £ = p3„, V ac , (which is between and 0.5). Further- 
more, it is a convex function of sin 2 £, which means that P ee is larger than 

the pure vacuum case for values of sin £ < - e % u vac . Away from the limit 

i+p ee ' 

P 6 L = cos 2 c<j, keeping in mind that the oscillatory terms average out, P ee is 
still larger than the pure vacuum case if cos 2 uj > sin 2 uj since P C L < cos 2 uj, 
as one can easily verify. 

Also, in the limit P^ = 1, sin 2 9 = 1/2, 



P =1 
w 2 



' Am 2 21L" 



2P, 



V 



;i - P c L )f/i + P C L U% + U% + 2^(1 - P C L )U^U, 2 cos 

V (4.17) 

The same expression applies for P TT with — > Z7 r ». This is a consequence of 
sin 2 6* = cos 2 6. Furthermore, in the limit sin 2 £ — > (and for sin 2 6 = cos 2 6), 
U^i = U Ti , which explains why P w = P TT for sin 2 £ < 10~ 2 . At VAC this 
equality remains for all values of sin 2 £. The reason for this is that, at VAC, 
the expression simplifies tremendously and P w = P TT = | (1 + cos 2 £). In 
the same region of the parameter space, the pure vacuum oscillation case 
yields P™ c = P™ c = |cos 4 £ — cos 2 £ + 1. Note that, in this region of the 
parameter space P™ c > P^, the inequality being saturated at cos 2 £ = 1. 

The same result also applies (approximately) at SMA, since the oscillatory 
terms are proportional to yj 'P C L (1 — P C L ) and 1 — P C L is very small at SMA (see 
Fig |3[). The equality P MM = P TT is broken at larger values of sin 2 £ because 
P C L 7^ cos 2 uj at SMA. 

It remains to discuss how P MAt and P TT diverge from the pure vacuum 
case at LMA and LOW. In the limit P C L = cos 2 a;, and averaging out the 
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oscillatory terms, 



Pw-P£ = ^ [ sin £ ( U % ~ (cos 2 ^^ + sin 2 ^ 2 2 )) - sin 2^ - U%) 

(4.18) 

This difference goes to zero as sin 2 £ — > 0. This is to be expected, since 
in this limit the difference of P 2 and cos 2 £ disappears. For small values of 
sin 2 £, the last term in Eq. f4.18| dominates, and, as discussed before, Uf tl — 
U* 2 = -O.5cos2cj + 0(sinO- Therefore, P^-P™ > (< 0) for cos2cj > 
(< 0). The expression for P TT can be obtained from Eq. ( |4. 18|) by replacing 

— > U T i and changing the sign of the last term. Therefore, since U^—U^ = 
-0.5 cos 2u + O(sin0, ^ - P£? > (< 0) for cos2cj < (> 0). When 
the oscillatory terms do not average out, it is easy to verify explicitly that 
the behaviour of the oscillatory terms follows the behaviour of the average 
terms, discussed above, and the inequalities obtained above still apply. 

The situation, however, changes, when P C L ^ cos 2 c<j, i.e., when matter 
effects due to the "M-L" system are relevant. In this region, a behaviour 
similar to the one observed in the "normal" hierarchy case is expected, since 
Am 2 ! > 0. Fig. [II] depicts constant P a p contours in the (Am^ x sin 2 em- 
plane. One should be able to see upon close inspection that the region 
P ee < 30% is smaller in Fig. [10] than the same region in the pure vacuum 
oscillation case, Fig |. Also, the constant P m (P TT ) contours are shifted 
to larger (smaller) values of sin 2 a;. The other prominent (and expected, as 
mentioned above) feature is the distortion of the contours at large values of 
Am 2 ^. This behaviour is similar to the one observed in Fig. 

I conclude this subsection with a comment on antineutrinos. As dis- 
cussed previously, F^^Am^Am^) = P a p{— Am| 1; — Am^), such that the 
"normal" hierarchies yield "inverted" hierarchy results for antineutrinos, and 
vice- verse. One cannot, however, apply Fig. ||| and Fig. [H] for the antineutri- 
nos because both Am 2 - have to change sign, not just Am 2 x . Qualitatively, 
however, it is possible to understand the constant P^q contours by examin- 
ing figures Fig. |8| and Fig. |TU| reflected in a mirror positioned at sin 2 u; = 0.5, 
meaning that P S( g(sin a; , Am 2 x ) ~ P Q/3 (cos 2 u, — Am^). The equality is not 
complete because one is also required to exchange 9 —>■ n — 6, as mentioned 
earlier. 
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Figure 10: Constant P a p contours in the (Am^ x sin 2 cj)-plane, at ATM (see 
text), for Am 2 ! < 0. 

4.4 Higher Neutrino Energies 

As the average neutrino energy increases, the values of P a p start to resemble 
more the pure vacuum case. This is easy to see from Figs. ^ and |3|. Any 
deviation of 1 — P 6 L from sin 2 u goes away even at LMA for E v ~ 50 GeV, 
while "H-M" effects remain important up to E v ~ 1 TeV, even though quan- 
titatively the effect decreases noticeably. This can be illustrated by the value 
of P% at ATM, for example, which drops from 0.87 for energies which range 
from 1 to 5 GeV (see the previous subsections) to 0.058, for energies which 
range from 100 to 110 GeV. 

Furthermore, all L° sc increase as the energy increases, for fixed values 
of Am 2 -. Therefore, LOW becomes indistinguishable from VAC at E v ~ 
100 GeV. For O(TeV) neutrinos the sensitivity to Am^ remains only for 
its highest allowed values, while one should start worrying about nontrivial 
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oscillatory effects due to L^ c . 

The case of higher energy neutrinos contains a more serious complica- 
tions: neutrino absorption inside the Sun. As the neutrino energy increases, 
one has to start worrying about the fact that absorptive neutrino interactions 
can take place. According to for neutrinos produced in the Sun's core, 
absorption becomes important for E v > 200 GeV. In this case, v e and 
interact with nuclear matter and produce electrons and muons, respectively. 
The former are capture and "lost" inside the Sun, while the latter stop before 
decaying into low energy neutrinos. The case of z/ r -Sun interactions is more 
interesting, because the r-leptons produced via charged current interactions 
decay before "stopping", yielding z/ r 's with slightly reduced energies. There- 
fore, it is possible to get a flux of very high energy initial state r-neutrinos 
but not muon or electron-type neutrinos. Such effects have been studied for 
high energy galactic neutrinos traversing the Earth 

The effect of neutrino oscillations inside the Sun in the presence of non- 
negligible neutrino absorption is certainly of great interest but is beyond the 
scope of this paper. 

5 Conclusions 

The oscillation probability of O(GeV) neutrinos of all flavours produced in 
the Sun's core has been computed, including matter effects, which are, in 
general, nontrivial. 

In particular, it was shown that, unlike the two-flavour oscillation case, in 
the three-flavour case the probability of a neutrino produced in the flavour 
eigenstate a to be detected as a flavour eigenstate (3 (P a p) is (in general) 
different from Pp a , even if the CP- violating phase of the MNS matrix van- 
ishes. This is, of course, expected since Sun-neutrino interactions explicitly 
break T-invariance. Indeed, it is the case of two-flavour oscillations which is 
special, in the sense that the number of independent oscillation probabilities 
is too small because of unitarity. 

The results of a particular scan of the parameter space are presented in 
Sec. 4. In this case, special attention was paid to the regions of the parameter 
space which are preferred by the current experimental situation. 

It turns out that, in the case of a "normal" neutrino mass hierarchy, it 
is possible to suppress P ee tremendously with respect to its pure vacuum 
oscillation values, by a mechanism that is similar to the well known MSW 
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effect in the case of two-ffavour oscillations: the parameters are such that 
electron-type neutrinos produced in the Sun's core exit the Sun (almost) 
as pure mass eigenstates, and the v e component of this eigenstate is small. 
Both P m and P TT can be significantly suppressed, and the constant P w 
and P TT contours as a function of the "solar" angle and the smaller mass- 
squared differences are nontrivial. One important feature is that when P w 
is significantly suppressed, P TT is not, and vice- versa. One consequence of 
this is that, for some regions of the parameter space, it is possible to have 
an enhancement of u T 's detected in the Earth with respect to the number of 
u^s (or vice- versa). This may have important implications for solar WIMP 
annihilation searches at neutrino telescopes, and will be studied in another 
oportunity. It is important to note that the effect of neutrino oscillations on 
the expected event rate at neutrino telescopes will depend on the expected 
production rate of individual neutrino species inside the Sun, which is, of 
course, model dependent. 

In the case of an "inverted" mass hierarchy, the situation is very similar to 
the pure vacuum case, and no particular suppression of any P aa is possible. 
Indeed, for a large region of the parameter space P ee is in fact enhanced, a 
feature which is also observed in the two-flavour case ||39|| . 

The case of higher energy neutrinos was very briefly discussed, and the 
crucial point is to note that, for neutrino energies above a few hundred GeV, 
the absorption of neutrinos by the Sun becomes important. The study of 
absorption effects is beyond the scope of this paper. 

Finally, it is important to reemphasise that the values of P a p computed 
here are to be understood as if they were evaluated at the Earth's surface. 
No Earth-matter effects have been included. It is possible that Earth-matter 
effects are important, especially the ones related to Am^, in the advent that 
{7g 3 = sin 2 £ turns out to be "large." 
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